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CORRELATION ANALYSIS AND SCINTILLATION FOR 
15-GHz LINE-OF-SIGHT PROPAGATION CHANNELS* 


by 

E. Mondre 

Goddard Space Flight Center 


INTRODUCTION 

Propagation studies of millimeter waves through the atmosphere are of concern both to com- 
munication applications and to investigations of the atmospheric structure. Most of the published 
results refer to long-term statistics of gaseous absorption and attenuation due to precipitation 
(References 1, 2, 3, and 4). There is a need for quantitative observation of short-term parameters 
such as correlation functions and the power spectral densities of the scintillations occurring on a 
line-of- sight path. These short-term statistical parameters allow a description of the time- and 
frequency- selective fading properties of a communication channel. In theoretical papers, dis- 
cussing error probability, interm odulation distortion (References 5 and 6) and optimal demodula- 
tion of analog signals (Reference 7), the assumption is made that channel characteristics can be 
accurately estimated. 

The millimeter-wave propagation experiment through the atmosphere between ATS-E space- 
craft and various ground stations throughout the United States will yield information about the feasi- 
bility of using such frequencies for communication links. The planned data reduction for this ex- 
periment includes long-term statistics of signal attenuation as well as channel characterization 
with correlation functions and spectral densities. The latter is referred to as short-term analysis . 
The amplitude of the received signal will not be a stationary variable if observed over a very long 
period (hours or days, for example). However, when observed over a short time (a few minutes) 
the observable quantities will, in general, be stationary and thus the system can be considered 
quasi- stationary. One has, in effect, a quasi-wide -sense stationary uncorrelated scattering chan- 
nel (References 8 and 9). It is the purpose of this report to discuss the short-term analysis based 
on experimental data obtained from a short line- of- sight link. 

EXPERIMENTAL ARRANGEMENT 

The path selected for the short-term analysis of a 15-GHz signal is over a distance of 4.6 km 
between a water tower in the Agricultural Research Center and Building 22 at Goadard Space 


•This work was accomplished while the author held a National Research Council Postdoctoral Resident Research Associateship sup- 
ported by National Aeronautics and Space Administration. 
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Flight Center. The length of the ground path is on the same order as that portion of a ground- 
spacecraft link which is influenced by atmospheric water vapor and oxygen as well as precipitation. 

The transmitting antenna is mounted at a height of approximately 30 meters on the water tower. 

The focal-point-fed paraboloid with 1.2-meter diameter has a nominal gain of 43 db and a half power 
beam width in the E -plane of 1.05 degree. Vertical polarization is employed throughout the experi- 
ment. The 15-GHz power source is a Varian Klystron with 760-milliwatt power output. 

On the receiving site at the roof of Building 22 (18 meters above ground) a horn antenna with 
an aperture of 40 cm 2 (nominal gain 19.5 db) is implemented. The signal is detected with a bal- 
anced mixer and a 30-MHz phase-lock loop receiver. The output of the receiver is recorded on a < 

strip chart recorder operating at a chart speed of 1 mm/sec and also on an FM magnetic tape re- 
corder. The profile of the ground along the path is fairly flat. Clearance of the line-of-sight path 
with respect to any natural obstacles is good. Note that in general the path and antenna alignments 
preclude ground reflection. 

Power and frequency instabilities of the klystrons are negligible since short-term analysis of 
atmospheric scintillations is the main goal of this experiment. During each of the short data runs 
(four minutes) power and frequency did not change. The receiver is calibrated on a clear day with 
low humidity with the aid of a calibrated attenuator. 

Rain data are collected at five points along the path. The rain gauges are spaced at approxi- 
mately equal distances (separation varies from 970 to 1240 meters). The resolution of the tipping 
buckets is 0.01 inch/tip. In determining the average rainfall rate along the path it is assumed that 
the indicated rate at a particular gauge applies over half the distance to the adjacent gauges. 


DATA RECORDING 

Data records are collected over a two-month interval (June and July). Selected samples for 
different atmospheric parameters (humidity, rain, wind) are used for data analysis. Two extreme 
cases are shown in Figure 1. Sample A is recorded on a clear day with few clouds, an average 
water-vapor concentration of 10 gr/m 3 , and a temperature of 70 C F. Record B is a sample taken 
during a rainstorm with a rainfall rate of 32 mm/hr, a water-vapor concentration of 17 gr/m 3 , 
and an average temperature of 76 C F. Record C is taken on a hot day (92°F) with a water-vapor con- 
centration of 16.5 gr/m 3 . The wind speed perpendicular to the path is approximately 11 km/hr for 
records A and B, and 20 km/hr for record C. 

The data on the analog magnetic tape are used in connection with a special-purpose multi- 
channel data processor built by Federal Scientific Corporation for the ATS-E Millimeter Wave 
Propagation Experiment. This processor can be used to determine correlation functions and spec- 
tra of the data records. (Details are discussed later.) However, the main part of the data analysis 
is performed on an IBM 360 digital computer. Therefore, the recorded data are played back through 
an analog-to-digital converter, to generate a digital data tape. The analog record is sampled at a 
30-Hz rate for 4 minutes. Each of the digital data records therefore contains 7200 samples. The 


2 



record- time length for each data runis selected 
to provide a reasonable trade-off between the 
statistically time-varying propagation medium 
and the statistical reliability of the measurement. 

The cutoff frequency due to sampling is 15 
Hz, which is well above all the expected fre- 
quency components of the scintillation. Ampli- 
0 tude quantization is done in 224 steps over the 

range of ±1 volt. This gives a resolution in 
¥ amplitude of about 9 mV. If a uniform probabil- 

ity distribution over one scale unit is assumed, 
the rms noise- to- signal ratio introduced by 
quantization is approximately 0.0013 (Refer- 
ence 10). 

STATISTICS OF SIGNAL FLUCTUATION 
Test for Stationary 

It is important to test the data records for stationarity, since calculation of correlation func- 
tion and spectrum is only useful for a channel which is at least a quasi-wide- sense stationary un- 
correlated scattering channel (References 8 and 9). To test the stationarity of the data, the sample 
record is divided into 20 equal time intervals. The data in each interval may now be considered 
independent. Both mean value and mean square value are computed for each interval. They are 
compared with the mean value and mean square value for the whole data record length. A run test 
(Reference 10) is used to decide whether the hypothesis (data are stationary) is accepted or re- 
jected. The test is performed at a 0.05 level of significance. 

During one occurrence with heavy rain the record length must be reduced to 1.5 minutes to 
ascertain stationarity of the data. All other records, even those with rain and strong wind, are 
stationary over the 4 -minute interval. 

This test verifies weak stationarity since mean value and mean square value only are tested, 
t But actual physical models are generally also strongly stationary if they are weakly stationary. In 

the case of data with a Gaussian probability density function, weakly stationary data are automat- 
ically strongly stationary. It is shown later that the expected Rician fading statistics can be ap- 
* proximated by a Gaussian distribution with nonzero mean. 

Signal Attenuation 

The mean signal level is calculated for each data run from 

* = (o 
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Figure 1— Samples of received signal amplitude. 
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where N is the number of samples, and x. are individual record samples. For data runs taken on 
humid days (water-vapor concentration between 10 and 30 gr/m 3 ) the mean value is approximately 
0.09 to 0.22 db/km below the dry-day level. This agrees with the reported values of 0.006 to 0.01 
db/km per gr/m 3 for attenuation due to water vapor (References 11 and 1). 

Data runs during precipitation are made in drizzle and rain with average rainfall rates of 16.5 
and 3 2 mm/hr. The readings of individual buckets during the heavy rain occurrence were not the 
same, thus indicating a spatial inhomogeneity along the path. The bucket readings varied from 2 to 20 
tips during a 5-minute period. The theoretical attenuation is calculated from the empirical relation 

A = a R b (2) 

where A is the attenuation in db/km, Rthe average rainfall rate in mm/hr, and a and b are param- 
eters. Numerical values for these two parameters can be found in the literature (References 3, 4, 
12, and 13). At a frequency of 15 GHz, suitable values seem to be a = 0.035 and b = 1.155. For 
the two reported rain occurrences we find attenuations of 0.9 and 1.92 db/km. The measured signal 
attenuation during these two events gives 0.9 and 2 db/km, which is in good agreement with the cal- 
culated values. This does not prove that the chosen values for parameters a and b will fit best for 
all rain events, since the attenuation depends in a very complex manner on the rainfall rate (drop- 
size distribution and rain intensity). The best choice for a and b can only be made from data col- 
lected over a long period during different rain events. But this is not within the scope of this paper. 

The calculated mean values are used to transform the data to have zero mean value. This is 
done in order to simplify subsequent formulas and improve the accuracy of the calculations. The 
data can now be converted into the ratio of received electric field to mean electric field by means 
of the receiver calibration curve. In the vicinity of the mean value the calibration curve is approx- 
imated by a straight line. The mean value m of the electric-field ratio E/E 0 is obviously equal to 
one. 

Signal Statistics 

The propagation channel can be described by theoretical models, a procedure that allows the 
determination of the expected probability density function for the received signal amplitude. The 
simplest model with the output of the channel a Gaussian random process leads to a Rayleigh dis- 
tribution of the envelope-detected signal. This model is usually not correct for line- of- sight prop- 
agation where a strong coherent (direct) signal is received. If the coherent signal is added to the 
Gaussian model, the envelope of the received signal will have a Rician statistic (Reference 14). 
Another model generally used in laser propagation studies assumes a logarithmic -normal (log- 
normal) distribution for the signal amplitude. This means that the logarithm of the received en- 
velope (or envelope expressed in db) is normally distributed. 

Rician and log-normal distributions are described by two parameters each, which can be cal- 
culated from the mean value and mean square value of the recorded data (Reference 15). Note that 
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the Rayleigh distribution is defined by one parameter only. It can be shown that in this case the 
ratio of mean value and standard deviation is constant and equal to 


m 

o 



1.92 . 


( 3 ) 


6 


r 

i 


All the data records have a ratio m/a greater than 25, which excludes, as expected, the Rayleigh 
distribution. For such large ratios of coherent and fading components, the Rician distribution 
can be approximated by a Gaussian distribution with nonzero mean (Reference 16). 


The probability density function for both models (Rician and log-normal) is computed for each 
data record. Figure 2 shows an example of the theoretical probability density for sample record B. 
It can be seen that there is a negligible difference between these two distributions, as is expected 
for such large ratios of steady and fading components. Figure 2 also shows the probability density 
calculated from record B, which is not unique because it is determined by the choice of the number 
of class intervals in which the whole range is divided (Reference 10). These are shown as "X"- 
marks in the figure. 


Figure 3 shows examples of the distribution 
function plotted on Gaussian probability paper. 
Again the Rician and log-normal distributions, 
as well as the experimental data, are plotted. 
The values at the "tails" of the distribution 



Figure 2 — Probability density function of received 
electric-field for record B. 



Figure 3— Rician, log-normal, and measured 
probability distribution. 
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function are statistically uncertain owing to the finite sample size. The "tail" values are derived 
from only a few data points; one additional data point falling into these intervals would cause a 
noticeable change in the distribution function. All data records can be classified as having either 
Rician or log-normal envelope density. 

Signal Scintillation 

Scintillation of the received signal is caused by fluctuation of the dielectric constant of the at- « 

mosphere. As already pointed out, ground reflections, which also would produce amplitude fading, 
are not to be expected. The standard deviation <j can be calculated from 



i= 1 


where x. are the transformed data values (electric-field ratios) with zero mean. For clear sky 
and calm wind the standard deviation of the field fluctuations ranges from 1 to 2 percent of the 
mean value, and goes up to3. 5 percent for heavy rain and wind velocities on the order of 20 km/hr. 

Since the amplitude is approximately a Gaussian distributed variable (see Figure 3), the prob- 
ability that the instantaneous signal level will be between +3~ and -3<? of its average level is 99.7 
percent. With this confidence level, the maximum scintillation of the signal amplitude in percent 
of the mean signal level is given by 


3ct 

A a = 4 — 100 , 


(5) 


where m = 1 and ^ is defined in Equation 4. The maximum scintillation is calculated for all recorded 
data runs. The absolute values for are between 3 and 10 percent. 

Tatar ski (Reference 17) developed a turbulence theory where the atmospheric turbulence is de- 
scribed by a "structure function." The structural constant in this structure function, denoted by C n , 
is related to the variance of the logarithm of the signal amplitude by 

a/ = 0. 31 C 2 k 7/6 L 11/6 

< n n 



for << i\T, « L 0 , where L is the path length and k r 2^/' the spatial wave number (\ = wavelength). 
L 0 is the outer scale of the turbulence (largest "blob" size). Close to ground it is often assumed 
that the outer scale is on the order of the beam height above ground. Near the ground, the inner 
scale i 0 is on the order of several millimeters. Equation (6) is derived on plane-wave assumption 
(Reference 17). 
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Although the condition i^Lis not exactly satisfied (L 0 ~ 20-25 m, Al = 9.6 m), values for 
C n are calculated from the measured signal fluctuation by means of Equation (6) and the definition 
of the variance of the logarithmic amplitude 


" ^ nE/E o) 2 ' ( 7 ) 

Since E/E 0 is approximately log-normal distributed, the quantity o-£ n can be derived from the al- 
ready computed variance of the signal amplitude a - 2 , according to 

= 'in (l + cr 2 /m 2 ) . (8) 

The structural constant C n is found to vary between 0.28 x 10” 6 and 0.93 x 10“ 6 m~ 1/3 . Reported 
values for C n obtained from measurements of fluctuations in phase differences and direct measure- 
ments to the refractive index and temperature spectra in the troposphere vary over a wide range 
from 0.019 x 10~ 6 to 0.98 x 10" 6 m" 1/3 (Reference 18). 

The maximum expected scintillation for a link through the entire atmosphere can be estimated 
if different models are assumed for the structural constant as a function of height. Tatarski T s 
theory gives the relation 




r c n 2 0 k"« 


( h o sec/3 ) 


11/6 


(9) 


t 

« 


where C n0 is the value of C n at the Earth's surface, and h Q is assumed to be 2000 m or 3200 m. The 
first value corresponds to the results of soundings with a balloon-borne refractometer (Reference 
19). The direction of the propagation link is defined by the zenith angle /3. The factor y depends 
on the assumed model. Values of y are shown in Table 1 (Reference 17). 


With Equations (5), (8), and (9) we obtain the expected maximum amplitude fluctuations sum- 
marized in Table 2. We calculate C n0 from the measured maximum value of C n = 0.93 x 10" 6 m" 1/; 
at an average beam height of 25 meters. A zenith angle of J3 = 48°30 T corresponds to the look angle 


from the ATS- E ground station at Rosman, North 
Carolina, to the ATS-E spacecraft in synchro- 
nous orbit at 108°W longitude. 


Table 1 

Values of y for Different Models for 
Structural Constant. 


Model 


7 

1 

C n 2 = C n 2 0 exp(-h/h 0 ) 

0.53 

2 

C n 2 = C/ 0 [l + (h/h 0 )3]- 1 

3.4 


Table 2 

Expected Scintillation for Earth -Spacecraft 
Link (15 GHz). 


Zenith 

Model 1 

Model 2 

Angle /3 

h Q = 2000 m 

3200 m 

2000 m 

3200 m 

0° 

±6.2% 

±9.6% 

±15.7% 

±24.3% 

48°30 1 

±9.1% 

±14.0% 

±23.0% 

±35.0% 


- 0.93 x 10"^ m at h - 25 m - 
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CHANNEL CORRELATION ANALYSIS 


The time- and frequency- selective fading properties of a communication channel can be char- 
acterized from the correlation function of a test signal (References 8 and 9). In this experiment 
only a single frequency is used; thus only the echo correlation function and its Fourier transform — 
the echo power spectrum— can be calculated. This leads to information about the coherence time 
and fading bandwidth of the radio channel. Two other parameters of interest are the coherence 
bandwidth and time spread of the channel. To measure those parameters we must transmit two 4 

carriers and evaluate the cross-correlation function as a function of their frequency separation. 

This is planned for the ATS-E Millimeter Wave Propagation Experiment. 

♦ 

Echo Correlation Function 

The echo correlation function is calculated from 


R(r) 


R( rh) 


1 

N- r 



x 


i + r 


r = 0, 1, 2 



( 10 ) 


where r is the lag number and h the sampling interval which is 1/30 second. The sample data x. 
are the transformed data with zero mean. It is very important to transform the data to zero mean 
before calculating the correlation function, in order to achieve accurate correlation function esti- 
mates. This is especially true if the mean value of the data is large compared with their standard 
deviation. For comparison it is better to use the covariance function rather than the correlation 
function. Dividing R(t) by R(0) gives the echo covariance function, which is shown in Figure 4 for 
different meteorological conditions. 


The echo correlation function is calculated from the envelope- detected signal rather than the 
complex signal amplitude. A relationship between envelope and complex covariance for Rician 



Figure 4— Echo covariance function of signal amplitude. 


envelope statistics is derived in Reference 15. 

As already pointed out, the ratio between coherent 
and random signal is very large. Therefore the 
calculated envelope covariance is indistinguish- 
ably close to the magnitude of the complex co- 
variance (see Reference 15). 

The coherence time T c of the channel is de- 
fined in the literature as the interval in time lag 
r over which R(r) is essentially nonzero. This 
definition of T c is useful from a theoretical point 
of view but it is usually not possible to calculate 
the correlation function with an accuracy for 
time lags great enough, to apply this definition. 


r 

♦ 
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For comparison and as a measure of T c (T c is used in an order -of- magnitude sense only) we ap- 
proximate R(t)/r(0) with the exponential function exp(-r/T c ), and redefine T c in terms of this re- 
lation. It can be seen from Figure 4 that T c is on the order of 1.5 to 4 seconds. For a physical 
interpretation of the coherence time one can think of T c as a measure of the average duration of 
fades on the channel. In transmitting digital data, the duration of a symbol should be considerably 
less than T c to allow for coherent integration of the received signal. 

Echo Power Spectral Density 

The power spectral density is obtained from the auto- correlation function by means of the 
Fourier transform. Owing to the finite sample size of the data records, only an estimate of the 
spectrum can be calculated. The auto- correlation function and power spectral density are, formally 
at least, even functions. Hence, the relation between those two functions may be expressed in the 
simpler form of the one-sided cosine transform, namely 

S(f) - 2 j R(r ) cos (2tt fr ) dr . (11) 

-'0 


A "raw" estimate of the true power spectral density is defined for an arbitrary f in the range 
0 < f <f c by 


S(f) 


2h 


R(0) 


n- i 

L 


R(rh) cos 


rrri 

T~ 


77n f 

R(nh ) cos ~f — 


(12) 


where f c = 1/2 h is the cutoff frequency and n <n/10. The spectrum is calculated for discrete fre- 
quencies f k = kf c / n ? with k = 0, 1, 2 • • • n. 

A final smooth estimate of the power spectral density may be obtained by frequency smoothing, 
which is performed by multiplying the auto-correlation function by suitable even functions of r m 
This amounts to weighting the correlation function nonuniformly. One way is to use the Hanning 
lag window weighting function, defined by 




D(rh) 


( 7TT 

1 + cos — 


l 0 


r = 0, 1 , * ■ * n 
r > n . 


(13) 


If time is discrete and computation digital, the smoothing can be done more economically in the 
frequency domain by transformation of the correlation function followed by convolution. Because 
D(rh) are finite sums of constants and cosines, so that their transforms are simple sums of Dirac 
delta functions with appropriate spacing, convolution means only multiplying with weights 0.25, 
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0.5, 0.25; viz. : 


s 0 = S(0) = 0.5S 0 + 0. 5 Sj , 

s k = s(k = 0.25S k _j + 0.5S k + 0.25S k+1 (k = 1 , 2, • • • n - 1 ) 

S n = S(f c ) = 0. 5 S n _ j + 0. 5 S n . 


( 14 ) 


Spectral distribution of the fluctuations indicate that for 15 GHz the scintillations are usually 
confined to frequencies less than 0.5 Hz. Figure 5 shows typical normalized echo power spectral 
densities obtained from the correlation functions of Figure 4. They are normalized with the square 
of the mean value. The ordinate in Figure 5 is the ratio of fluctuation power density and mean power. 
Plotted are the two extreme conditions A and B (see Figure 1) with standard deviations of 0.01 and 



Figure 5 — Echo power spectral density of amplitude 
fluctuations (normalized with mean d.c. power). 


0.034. Also plotted is the spectrum of record C 
with a standard deviation of 0.025. 

Since the total fluctuation power P ac is equal 
to the variance o 2 9 that is 

p p ^2 

S(f)df = pf = (15 ) 

J 0 dc m 

(S(f ) is normalized with square of the mean 
value), we may normalize the power spectral 
density by multiplying the plots in Figure 5 by 
m 2 /a 2 . This leads to a presentation where the 
area under each curve is equal to one. As shown 
in Figure 6, all data records show a very similar 
power spectral density. The ordinate is now the 
ratio of fluctuation power density to total fluctua- 
tion power. 

The fading bandwidth of the 4.6 km communi- 
cation channel is approximately 0.4 Hz and does 
not seem to change very much with weather 
conditions. 

It is interesting to compare the power spec- 
tral density with the theoretically expected spec- 
trum derived from the spatial correlation and at- 
mospheric structural functions. According to 
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Reference 18 the spectrum of the amplitudinal 
fluctuations is described by 

S(f) = 0.016 C n 2 k 2 Lf 8/3 v T 5/3 , (16) 

where v T is the component of the mean wind 
speed normal to the path. This equation is valid 
for frequencies f > v T /V\L whereas for f < v r /Ykh 
the spectrum is expected to be nearly constant. 

Figures 7 and 8 show examples of the 
measured as well as the theoretical spectrum 
for data records A and C. In Figure 7 there is 
good agreement for f/v T less than 0.5 m" 1 (f 
= 1.5 Hz). The deviation for f/v T > 0.5 is due to 
noise. The fluctuation power is concentrated in 



f/v T (m' 1 ) 

Figure 7— Normalized spectrum of amplitude 
fluctuation for record A. 



FREQUENCY (Hz) 

Figure 6— Echo power spectral density of amplitude 
fluctuations (normalized with total fluctiation power). 



nr 3 nr 2 io _l io° io 1 


f/v T (m-1) 

Figure 8— Normalized spectrum of amplitude 
fluctuation for record C. 


frequencies below 0.5 Hz; therefore for higher frequencies the noise contribution predominates. 
Assuming white noise, the theoretical curve in Figure 7 should have a +45-degree slope for higher 
frequencies. The fluctuation power for record C is well above noise level for frequencies up to 
6 Hz, as can be seen in Figure 8. 
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Analog Correlation Analysis 


For comparison, the channel correlation analysis is also carried out on a special-purpose 
multi-channel data processor built by Federal Scientific Corporation for the ATS-E Millimeter 
Wave Propagation Experiment. A description of this processor is given in Reference 20. As in- 
put for the data processor, the analog tape with the recorded receiver output voltage is used. 

The auto- correlation functions are calculated from data records 153.6 seconds in length. The 
time lag of the correlation function is from 0 to 102.4 seconds, which is up to 2/3 of the record 
length. However, only small time lags are of interest, since the coherence time T c is on the order 
of a few seconds. The data processor can perform the correlation analysis from shorter record 
lengths. However, in order to minimize the statistical uncertainty of the correlation function, the 
record length should be much greater than T c . The obtained values of T c for the investigated data 
records are comparable with those already discussed in the digital analysis. 



Figure 9— Power spectral density vs. frequency for 
records B and C. 


The spectrum calculation is performed over 
the 15-Hz analysis range. Again, all records are 
analyzed and compared with the results of the 
digital data reduction. Examples for data records 
B and C are shown in Figure 9. 

The record length used for computation in 
the 15-Hz range is 30 seconds. The resolution 
bandwidth of the spectrum is 0.0325 Hz. Since 
all the power spectral densities are derived 
from data samples, they are statistical esti- 
mates of the true spectrum. The latter could 
be obtained by using certain smoothing and aver- 
aging techniques. The ordinate axis in Figure 9 
is not calibrated. The derived fading bandwidths 
of the data records are approximately equal to 
those obtained in the digital analysis. 

The peak at 3.7 Hz, which is present in all 
of the analog and digital derived power spectral 
densities, is believed to originate from mechan- 
ical resonances of the transmitter antenna. 


CONCLUSIONS 

The outlined data reduction and experimental data for the short line- of- sight link are a basis 
for comparison with the data obtained from the ATS-E Millimeter Wave Propagation Experiment. 
It has been shown that a Rician fading channel is a proper model for a line- of- sight path. 
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Numerical values for the channel parameters are obtained. The coherence time is on the 
order of 1.5-4 seconds, and the fading bandwidth is below 0.4 Hz. The main influence of precipita- 
tion on wave propagation at 15 GHz is attenuation. No noticeable change of the fading bandwidth 
due to rain was observed. 

It is important that correlation analysis should be performed when the variance of the signal 
is above a certain threshold (determined by the sensitivity and accuracy of the experiment). As 
expected, fairly large scintillations of ±1 db are observed during increased windspeeds. The ob- 
tained channel parameters (coherence time, fading bandwidth, etc.) could be used for investigating 
long-term changes. To achieve this, the experiment must be performed over a much longer period 
(one year). 

The estimated maximum scintillation for a link throughout the atmosphere depends on the 
chosen theoretical model, and can be verified only by an actual measurement. The obtained spec- 
tral density of the scintillations for the line-of-sight link are in good agreement with the theoreti- 
cally derived fluctuation spectra. 
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